Unbalanced instabilities of rapidly rotating stratified shear flows 



J. Vanneste 1 and I. Yavneh 2 

1 School of Mathematics and Maxwell Institute for Mathematical Sciences, University of Edinburgh, 

Edinburgh EH9 3JZ, UK 

2 Department of Computer Science, Technion, Haifa 32000, Israel 



Abstract 

The linear stability of a rotating, stratified, inviscid horizontal plane Couette flow in a channel is 
studied in the limit of strong rotation and stratification. Two dimensionless parameters characterize 
the flow: the Rossby number e, defined as the ratio of the shear to the Coriolis frequency and assumed 
small, and the ratio s of the Coriolis frequency to the buoyancy frequency, assumed to satisfy s < 1. 
An energy argument is used to show that unstable perturbations must have large, 0(e _1 ) wavenumbers. 
This motivates the use of a WKB-approach which, in the first instance, provides an approximation for 
the dispersion relation of the various waves that can propagate in the flow. These are Kelvin waves, 
trapped near the channel walls, and inertia-gravity waves with or without turning points. 

Although, the wave phase speeds are found to be real to all algebraic orders in e, we establish that the 
flow is unconditionally unstable. This is the result of linear resonances between waves with oppositely 
signed wave momenta. Three modes of instabilities are identified, corresponding to the resonance between 
(i) a pair of Kelvin waves, (ii) a Kelvin wave and an inertia-gravity wave, and (iii) a pair of inertia-gravity 
waves. Whilst all three modes of instability are active when the Couette flow is anticyclonic, mode (iii) 
is the only possible instability mechanism when the flow is cyclonic. 

We derive asymptotic estimates for the instability growth rates. These are exponentially small in e, 
of the form lm lu — aexp( — \&/e) for some positive constants a and For the Kelvin- wave instabilities 
(i), we obtain analytic expressions for a and ty; the maximum growth rate, in particular, corresponds 
to \l/ = 2. For the other types of instabilities, we make the simplifying assumption s C 1 and find that 
\& = 2.80 for (ii) and $ = n for (iii). The asymptotic results are confirmed by numerical computations. 
These reveal, in particular, that the instabilities (iii) have much smaller growth rates in cyclonic flows 
than in anticyclonic flows, in spite of ha ving both >fr = tt. 

Our results, which extend those of iKushner et all (| 19981 ) and lYavneh et ah! l|200ll ). highlight the 
limitations of the so-called balanced models, widely used in geophysical fluid dynamics, which filter out 
Kelvin and inertia- gravity waves and hence predict the stability of the Couette flow. They are also 
relevant to the stability of Taylor-Couette flows and of astrophysical accretion discs. 



1 Introduction 

Rapid rotation and strong density stratification characterise the dynamics of geophysical fluids, the atmo- 
sphere and the oceans in particular. Two dimensionless numbers are used to measure the importance of 
these two effects relative to nonlinear advection: the Rossby number 

U 

and the Froude number 

f = JL. 

ND 

Here U is a typical horizontal velocity, / > is the Coriolis parameter, N the Brunt -Vaisala frequency, and 
L and D are typical horizontal and vertical length scales. With N > /, as is realistically the case, the Rossby 
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Figure 1: Schematic of the Couette flow, with velocity (Ay, 0, 0), of a three-dimensional Bousinesq fluid with 
constant Brunt -Vaisala frequency N. The domain rotates around the (vertical) z-axis at rate //2; it is 
unbounded in the x and z directions, and bounded in the y directions by two rigid walls at a distance 2L. 



number estimates the maximum ratio between the typical frequency of the (slow) advective motion (given by 
U/L), and the frequency of inertia-gravity waves (bounded from below by /). Its smallness, explicity e«l, 
has an important dynamical consequence, namely the weakness of the interaction between advective motion 
and inertia-gravity waves. This, together with the observation that inertia-gravity waves have generally 
weak amplitudes in the atmosphere and oceans, has led to development — and success — of the so-called 
balanced models, which filter out inertia-gravity waves completely. These models describe only the slow, 
large-scale dynamics, termed balanced because of its closeness to hydrostatic and geostrophic balance. They 
can be derived asymptotically, using power-series expansions in e, and in principle can achieve an arbitrary 



can be derived asymptotically, using power-series expansions 
algebraic accuracy 0(e n ) re.g.. lWardfl997lrWarn et al.lll995l) 



To understand balanced dynamics and its limitations more fully, it is important to identify and quantify 
the phenomena that balanced models fail to capture. Of particular interest are those unbalanced phenomena 
which occur in spite of the smallness of e and cannot be suppressed by balancing the initial data. In the 
present paper we consider one such mechanism, namely the instability of balanced flows to unbalanced, 
gravity-wave-like perturbations. Since this type of instability is absent from balanced models of arbitrary 
high accuracy (which all have qualitively similar stability conditions; see lRen &: Shepherd 1997 ). the growth 



rates can be expected to be o(e n ) for all n > 1 or, in other words, to be beyond all orders in e, and typically 
exponentially small in e. Our results confirm this scaling and show that the instability bands, i.e., the range 
of unstable wavenumbers, are exponentially narrow. 

We note that unbalanced instabilities like the one examine d in this paper are di s tinct from the mechanism 
of spontaneous generation of inertia-gravity waves sudied in IVanneste fc Yavneh ( 2004 ). Both mechanisms 



are exponentially weak, but whilst the exponentially small quantity is the growth rate for instabilities, it 
is the amplitude of the waves in the case of spontaneous generation. This difference may not be essential, 
however, if the unbalanced instabilities saturate at a level that decreases to zero with growth rate, as is 
typical. Another difference is the fact that the instabilities require an initial unbalanced perturbation, 
whilst spontaneous generation occurs from entirely balanced initial conditions. We emphasize that both 
mechanisms provide potential sources of inertia-gravity waves in the atmosphere and oceans. What the 
exponential smallness indicates in both cases is that the effectiveness of these sources is highly sensitive to 
the Rossby number. 

The specific flow whose stability we study is a horizontal Couette flow with velocity (Ay, 0,0), modelled 
using the Boussincsq approximation with constant N, and an /-channel of width 2L. See Figure [T] for an 
illustration. A natural definition of a (signed) Rossby number for this flow is the ratio 

A 

e = 7 
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of (minus) the basic-flow vorticity to the planetary vorticity. For e > (< 0), the shear is anticyclonic 
(cyclonic). The other dimensionless parameter characterising the flow can be taken to be the Prandlt ratio 

/ 

We restrict our attention to s < 1 and note that in the atmosphere and oceans s«l generally holds. 

Because it is steady, the flow under consideration remains exactly balanced for all times, unlike generic 
time-dependent flows. Furthermore, it is stable in any balanced approximation, however accurate: this 
is because the shear is linear, and hence the potential vorticity constant, whilst balanced instabilities are 
inflectional instabilities, which require changes in the sign of the potential-vorticity gradient. Thus, with this 
flow, there are none of the difficulties in separating inertia-gravity waves from balanced motion that would 
appear for more complicated flows, and the analysis reduces to a straighforward linear stability analysis. 
The smallness of e is of course exploited to derive asymptotic results. 

A number of authors have investigated gravity-wave-like instabilities of shear flows, although mostly in 



geometry (Satomurab 


.98ldd 


Dritschel & Vannestd 


2006), 



The emphasis was not, however, on the small e l imit; indeed , in shallow water, flows with e<Cl and F = 0(e) 
are linearly stable as Ripa's theorem indicates (|Ripa Il983h . In contrast, the three-dimensional model exam- 
ined here turns out to be always unstable, with growing modes whose horizontal and vertical wavenumbers 
scale like e _1 . Our analysis has nevertheless many common features with some of the works cited above, in 
particular the use of the WKB approximation. A common theme (in particular with Naravan et al.lb~987t ) 
is also the interpretation of the instabilities in terms of (linear) resonances between modes with differently 
signs of the conserved w ave energy (or pseudoenergy) and wave momentum (or pseudomomentum) (see, e.g., 
Craik Il985l iRipal 1990l and references therein) . 

In the presence of lateral boundaries, as is the case here, there are two types of unbalanced modes: inertia- 
gravity waves, which are oscillatory in the cross-stream direction, and Kelvin waves, which are trapped at each 



bound ary. Instabilities involving the resonance of Kelvin waves h ave been studied recently bylKushner et al 



( 19981 ) for the model considered here, and bv lYavneh et al. ( 200lh and Molemaker et akl (12001) in the annular 
geometry of the (stratified) Taylor-Couette flow (see also iRiidiger et al. I (|2002h and iDubrulle et all (2005) 
for astrophysica l applications) . For simpl e geometric reasons, the se instabilities occur only for anticyclonic 
shears (A > 0). lYavneh et all (|200ll ) and iMolemaker et all (j200lh also identified other modes of instability 
in anticyclonic shears. These can be associated with the resonance between Kelvin and inertia-gravity 
waves, and b etween inerti a -gravity waves. The first mechanism is analogous t o the mixed-mode instabilities 
examined bv lSakail (|l989h . iMcWilliams et all (|2004 ) . IMolemaker et all (|2005h and IPlougonven et all (|2005l ) 
in a variety of contexts. As we show, the second mechanism is also active in cyclonic shears (A < 0). Thus, 
we establish that the stratified horizontal Couette flow is unconditionally unstable. 

For all the instabilities that we study, the growth rates are exponentially small in e because the resonant 
waves with differently signed wave momentum are localised exponentially in different sides of the channel. 
We provide both a qualitative description of the instabilities, based on the mode resonance and conservation 
laws, and quantitative results, based on the WKB approximation and numerical computations. 

The remainder of this paper is organized as follows. The linearized equations of motion governing the 
evolution of perturbations in the Couette flow are introduced in Sj2l The conservation laws for the wave 
momentum and wave energy are also introduced there. The latter conservation law is used to show that the 
horizontal and vertical wavenumbers of growing perturbations must be 0(e~ 1 ) or larger. This motivates the 
WKB approach developed in §jj3]-|l| In Sj3]we formulate the eigenvalue problem for the normal modes of the 
system, then provide an approximate solution using a WKB expansion f £|3.ip . To all orders in e, this leads to 
purely real eigenfrequencies or, in other words, to waves rather than growing modes. Instabilities with growth 
rates beyond all orders in e are however possible, and we go on to show that they do occur. Focusing on the 
modes susceptible to be involved in instabilities, we give some details of the dispersion relation and structure 
of Kelvin waves ( £|3.2p and inertia-gravity- waves with turning points f £|3.3p . We then use arguments based 
on wave-momentum signature to show that the linear resonance between waves does lead to instabilities for 
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both cyclonic and anticyclonic shears f H3.5j) . Section Q] is devoted to the estimation of the instability g rowth 
rate s. A detailed a symptotic estimate for Kelvin- wave instabilities, extending those of lKushner et al. I (|l998l) 
and lYavneh et aT. (2001), is derived in £)4.1I Rough estimates (focusing on the exponential dependence and 
ignoring order-one prefactors) are then obtained for the weaker types of instabilities (§ i)4.2tfO]) . These 
estimates are confirmed by the numerical solutions of the eigenvalue problem presented in £14.41 The paper 
concludes with a Discussion in SJS1 



2 Model 



We consider small-amplitude perturbations to the Couette flow described in the Introduction and in Figure 
[TJ The corresponding linearized equations of motion can be written as 



D t u-(f-A)v 


= -d x p, 


(2.1) 


D t v + fu 


= -dyP, 


(2.2) 


D t w + p 


= -d z p, 


(2.3) 


Dtp- N 2 w 


= o, 


(2.4) 


d x u + d y v + d z w 


= o, 


(2.5) 



where (u, v, w) are the components of the velocity perturbation, p is pressure perturbation, p the buoyancy 
perturbation, and D t — dt + Ayd x . The material conservation 



D t q = 



of the perturbation potential vorticity 



q= (/- K)d zP -N 2 {d x v-d y u) 

follows readily. We restrict our attention to pert urbations with vanishing po tential vorticity, q = 0, since 
this is a characteristic of unbalanced motion. (See IVanneste fc Yavneh ( 2004 ) for a study of the generation 
of inertia-gravity waves from perturbations with q ^ 0.) With this restriction, the conservations of the wave 
energy (pseudoenergy) 



£ 



|u£ 

2 



2N 2 



Ay 



ud z p - wd x p 



N 2 



dxdydz 



and of the wave momentum (pseudomomentum) 



M = 



ud z p - wd x p 
N 2 



dxdydz 



(2.6) 



(2.7) 



are readily derived, as detailed in Appendix lAl 

The conservation of £ constrains the structu re of unstab le perturbations. This is because exponentially 
growing modes must have vanishing £ (see, e.g.. iRipa 1990f) . Completing the squares in (|2.6|) , we rewrite £ 
as 



« - \ 



Ayd z p 
N 2 



Ayd x p 
N 2 



2„,2 



pr A 2 y . 2 



{dzpf 



dxdydz. 



Clearly, instability can only occur if the perturbation satisfies A 2 y 2 \d x p) 2 + (d z p) 2 ] > p 2 N 2 somewhere in 
the channel. In terms of horizontal and vertical wavenumbers k and m, this gives the condition 



N 



(k 2 



2)1/2 



< AL, 



(2.8) 
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which can be recognized as a subsonic condition: instability occurs only for modes whose phase speed is less 
than the maximum basic- flow velocity. With s < 1 as assumed, the subsonic condition implies that L{k 2 + 
m 2 ) 1//2 > e _1 , and therefore that modes involved in instabilies have asymptotically large wavenumbers. One 
interpretation of this result states that the Rossby number based on the wave scale, that is, h.L{k 2 +m 2 ) 1 / 2 / f , 
is greater than unity for unstable modes. 

We note that for the shallow-water model with d epth H, the subsonic condition analogous to (|2.8p is 
(gH) 1 / 2 < AL and does not involve the wavenumbers ( Ripa 19831 ). It is never satisfied for sufficiently small 
A and thus, for order-one Burger number, for sufficiently small e. Thus the shallow- water analogue of our 
model is linearly stable in the limit e — ► 0. 



3 Normal modes 

Let us now consider normal- mode solutions of the linearized equations of motion (|2.ip - (|2.5p . The subsonic 
condition (12. 8|) suggests that the wavenumbers k and to should be rescaled by e. We therefore write the 
dependent variables in the form 

u(x, y, z, t) — u(y/L) exp [ie~ 1 L~ 1 (kx + mz/s) — fuot\ , (3-1) 

with similar expressions for v,w,p and p. Here k, m and to are dimensionless wavenumbers and frequency, 
with their dimensional counterparts given by k/(eL), m/(esL) and feu, respectively. Without loss of gen- 
erality we assume that k > 0. Note that the non-dimensionalisation then implies that modes with lu > 
(lu < 0) propagate to the right (left) in anticyclonic shear and to the left (right) in cyclonic shear. 
In terms of the dimensionless k and to, the subsonic condition (|2.8p reads 

r = (s 2 k 2 + m 2 ) 1 / 2 > 1. (3.2) 

Introducing the normal modes (|3.ip into (|2.ip - (|2.5p leads to a system of ordinary differential equations 
for u, v, w, p and p. These independent variables can be eliminated in favour of p, leading in particular to 

1 e(l — e)p' + kujp „ lN m „ 1 muj „ 

jP and w— — — jnrP: (3-3) 



efL u; 2 -l + e ' r e/Ll-s 2 ^ 2 " eNL 1 

where prime denotes differentiation with respect to the dimensionless variabl e y/L which we henc eforth 
denote simply by y. A second-order differential equation, already obtained by iKushner et al. then 
follows. It reads 

e¥' - j^^P* - ( e l -^±± + to 2 ±-^ P = 0, (3.4) 
where 

uj = lu — ky. 

It is supplemented by the boundary conditions v — 0, that is, 

ecp' +p = at y = ±1, (3.5) 

where c = c — y = oj/k — y. Note that the singularities of (|3.4p for Co 2 = 1 — e are removable: in particular, 
they are absent from the equation for u equivalent to (|3.4p and given in Appendix [Bj 

3.1 WKB approximation 

Together, (|3.4p and (|3.5p constitute an eigenvalue problem from which the dispersion relation giving w as a 
function of k and to can be derived. Taking advantage of the small parameter e, this eigenvalue problem can 
be solved approximately using the WKB method. To this end, we first expand (|3.4p in powers of e, with the 
frequency 

lu = ujq + euji + • • • 
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turning out to be real to all orders. Taking into account that p' = 0(e 1 ) and p" = 0(e 2 ), we rewrite (|3 .4[) 
as 

2e 2 ku} ~, 



2 ~// 
e p 



\ 2 p- 



p' + ehp = 0{e*), 



where 



and 



A = k 2 



2 1-^0 

1 



(3.6) 
(3.7) 



2fc 2 



1 1 - S 2 CJ 2 



We introduce solutions of the form 

p = (g± + m± 

and find that g± satisfies 



) exp 



■ 2wi 



m 2 (l — s 2 )o>o 



(1 - s 2 ^ 2 ) 2 



Kv') <V 



into 



4 
.9± 



A' 



kLu Q 



2A ' 1 ' 



7i 

2A' 



(3.8) 



(3.9) 



where cr = sgne equals +1 for an anticyclonic shear and —1 for a cyclonic shear. The solution can be written 
as 



9± = A 



A 



2\ 1/2 



exp 



To - 



2A(y') 



(3.10) 



where A is an arbitrary complex constant. Note that this solution is single- valued near ljq = ±1, consistent 
with the observation that the singularites of (|3.4[) for Co 2 = 1 — e are removable: the multi-valuedness caused 
by the square root factor in (|3.10[) is cancelled by that of the integral in the argument of the exponential. 
We can classify the solutions (|3.8| according to the sign of A 2 in the channel and distinguish: 



Kelvin waves (KWs) , for which A 2 > for — 1 < y < 1. These modes are trapped exponentially near 
one of the boundary, with 0(e) trapping scale. 

Inertia-gravity waves (IGWs) , which satisfy A 2 < in at least part of the channel. There they have an 
oscillatory structure with O(e) wavelength. 

We now derive approximate dispersion relations for both types of waves. Together with information on 
the signature of their wave momentum discussed in £13.51 these allow the prediction of instabilities associated 
with KW-KW, KW-IGW and IGW-IGW resonances. Asymptotic estimates for the growth rates of these 
instabilities are derived in 21 where they are compared with numerical results. 



3.2 Kelvin waves 

We first consider WKB solutions to f|3.4|) for which A 2 > 0. Two independent solutions can be written as 



P- 



1-Cof 



2\ 1/2 



exp 



2 \ I/ 2 



exp 





Kv') - 


%') ' 

£ 2A(y') 


dy' 


k 


fO(6)] 


-7; 


Kv') - 


%') 1 
2A(y')J 




>[n 


-0(e)] 



The dispersion relation is found from the boundary conditions in the form 

ec(-l)p'_(-l) +p_(-l) ec(-iy + (-l) +p+(-l) 
ec(iy_(l) +p_(l) ec(iy + (l) +p+(l) 



0. 



(3.11) 
(3.12) 

(3.13) 



G 



Since the off-diagonal terms are exponentially small, the dispersion relation factorises to all orders into two 
branches corresponding to KWs trapped at each boundary. We denote by KW± the branch trapped at 
y = ±1, respectively; the corresponding frequency satisfies 

ec(l)p' + (l) +p+(l) = and ec{-l)p'_{-l) = 0. (3.14) 

At leading order in e, these two relations reduce to 

ctc (1)A(1) + 1 = and - crco(-l)A(-l) + 1 = 0. 

with solutions 

c (l) = -- and co(-l) = -. (3.15) 
r r 

(In addition, there are spurious solutions kco — ±o\) Thus, the KWs localised near y = ±1 have the 
leading-order dispersion relation 

Co = 1 and c = — 1 H — . (3.16) 

r r 

Higher-order approximations for the KW dispersion relation can be obtained by pursuing the expansion in 
powers of e, each leading to a purely real correction to (|3.15|) . 



3.3 Inertia-gravity waves 

In the region where the IGW is oscillatory, two independent WKB solutions (|3.8| can be written as 

1/2 



P = 



1-CM 



cxp i ±i\e\ 



e w) dy . 



(3.17) 



where t > is defined by 



I 2 = -A 2 . 



Depending on the value of w, IGWs can have at most two turning points, i.e. points where A = £ — 0, in the 
channel. These are located at 



y± = cq 



±-(l + 8 2 ) 1 / 2 , where 



5 = m/k, 



(3.18) 



on either side of the 'critical level' y = Co where luq = 0. The mode structure is then oscillatory for y < y_ 
and y > y+, and exponential for ?/_ < y < y+. Here, we concentrate on modes with at least one turning 
point since, as argued in ^3.51 below, the presence of a turning point is necessary for instability. These IGWs 
are localised on one side of the channel and exponentially small on the opposite boundary. 

Let us consider one such IGW that is decaying exponentially with y in [y-,y+] and denote the corre- 
sponding solution by (Its counterpart, growing exponentially in [y_,y + ] and denoted by p+, is readily 
deduced using the symmetry (y, c) i— > (— y, — c).) In [y_,y + ], the solution p_ can be written as 



A 



1-^0 
A 



1/2 



exp ■ 



A(y')-e 



%') 



2A(y') 



dy' 



(3.19) 



The boundary condition (|3.5[) at y = 1 is satisfied automatically to all orders in e. The form (|3 . 1 9[) breaks 
down in an e 2 / 3 neighbourhood of where it is replaced by an Airy function Ai. In [—1, y_], the solution 
is given by a linear combination of the two solutions in (|3.17l) . The connection formula, which relates the 
two arbitrary constants to A and is found by matching with the Airy function, gives (cf. Bender & Orszag, 
Eq. (10.4.16)) 

p _~2a( 1 —^) sinjler 1 ^" i{y') 



I 



%') 



2l{y<) 



Ay' 



(3.20) 
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The dispersion relation is then found by applying the boundary condition (13.14[) at y = — 1, leading to 



where 



Solving for S(y-), we find 



crc(- l)£(-l) cos %_) + sin 



0(e), 



5(2/-) 



'2*(y') 



dy' 



7T 

4' 



5(y_) = mr + tan- 1 [crc(-l)^(-l)] + 0(e) 



where n is an integer. At leading order this gives 

rv- 



£(y) dy 



nn e 



(3.21) 



which determines Co implicitly. The next order relation determines c±. 

Let us write the dispersion relation (|3.21|) for cq in a convenient form. Define \x by 



Co 



-1 



(l + J 2 ) 1 ^ 



(A* + 1) 



where 5 — m/k, so that 



y- 



-l 



(l + J 2 ) 1 ^ 



-/'• 



The assumption that this turning point is inside the channel imposes the restriction < /i < 2r(l + S 2 ) 1 ' 2 . 
Introducing the integration variable F, with y = — 1 + (1 + 5 2 ) 1 / 2 y/r, reduces Q3.2ip to the expression 



1 



i -i/ 2 (y-^i - l) 2 



1/2 



dr = 



(s 2 + 5 2 ) 1/2 
1 + S 2 ' 



(3.22) 



where 



s 2 {l + S 2 ) 



s 2 + S 2 

This defines implicitly a function s, n|e|) with values in [0, v^ 1 ~ 1], from which Co is deduced. Taking 
both the solution p_ and its symmetric p+ into account, we find the two branches 



c = ±1 =F 



(1 + S 2 ) 1 / 2 



H<5,s,n|e|) + 1], 



(3.23) 



corresponding to modes exponentially small near y = =pl and denoted by IGW±, respectively. Again, higher- 
order approximations to the phase velocity can in principle be computed, leading to real corrections to cq in 
powers of e. Note that, at leading order in e, the dispersion relation is the same for both signs of e, that is, 
for both cyclonic and anticyclonic flows. An asymmetry only appears at higher order. 
For n — 0(1), \i — > as e — + 0, and the leading-order dispersion relation reduces to 



Co 



±1T 



(1 + 5 2 ) 1 ' 2 



(3.24) 



corresponding to y± — > ±1. The small- /i behaviour of the left-hand side of (|3.22p then suggests that the 
successive branches n = 1, 2, • • ■ are 0(e 2 / 3 ) apart. 

The asymptotic results (|3 . and (|3.23|) provide a first approximation to the dispersion relation of KWs 
and IGWs. We have extended this by solving the eigenvalue problem (|3.4[) — (|3.5[) (or rather the equivalent 
formulation (IB. II) - (IB. 2]i in terms of u) numerically Our numerical solver is the same as the one used in 



ED 

D d20C 



Yavneh et al.l (|200lh . employing a second-order finite- volume discretization of (|B.1|) - (|B.2[) . For given physical 



parameters and wavenumbers, m and k, we search for eigenfrequencies for which the matrix representing the 
discretized system is singular. The codes are implemented in MATLAB, with the search performed using 
the fminsearch function that employs the so-called Simplex algorithm. 
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Figure 2: Dispersion relation for an anticyclonic flow, with e = 0.1, 5 = 2 and s = 0.1. The two Kelvin waves 
(labelled KW) are shown along with many inertia-gravity waves (labelled IGW) . The asymptotic estimates 
for |e| <C 1 (solid curves) are compared with numerical solutions (dotted curves). 



3.4 Dispersion relation 

Figures [2] and [3] show the dispersion relation for anticyclonic and cyclonic flows, respectively. The parameters 
have been chosen as e = 0.1, 5 = 2 and s = 0.1, but the qualitative features remain the same for a wide range 
of values. The numerical results (dotted curves) are compared with the asymptotic estimates (solid curves) 
to confirm the validity of the latter. For KWs, we have used an 0(e)-accurate estimate which improves on 
(I3.16[) by adding the term eci = =Feer(2r) derived in Appendix [Cl For IGWs, we have used the estimate 
(I3.23[) . corrected in the anticyclonic case by subtracting e/2 from the square bracket. This correction, which 
can be viewed as an experimentally determined O(e) term in the expansion of c, is made for the clarity of 
the plot: without it, the O(e) error in the dispersion relation is not significantly smaller than the 0(e 2 / 3 ) 
distance between branches, and it is difficult to relate each asymptotic curve to its numerical counterpart. 
No corrections were necessary for the cyclonic shear, suggesting the dispersion relation (|3.23p is already 
0(e)-accurate in this case. To confirm this would require to continue the asymptotic developments to the 
next order in e; this is a daunting task which we have not attempted for IGWs. 

Figures [2] and [3] demonstrate the multiple intersections between the branches IGW± of the dispersion 
relation. In the anticyclonic case, there are additional intersections between the branches KW±, and between 
KW± and IGW T . (The KW do not appear in Figure[3]for the cyclonic case because they have |c| > 1.) The 
intersections, associated with the linear resonance between modes, are generically spurious: they result from 
the finite resolution of the plot for the numerical results, and from the limited accuracy for the asymptotic 
ones. There are in fact two possible behaviours: (i) mode conversion, when the phase velocities remaining real 
and the two curves, rather than intersecting, locally form the two branches of a hyperbola, or (ii) instability, 
when the phase velocities on the two branches become complex conjugate with non-zero imaginary parts. 
The two situations are distinguished by the signs of quadratic invariants, such as the wave momentum, along 
the colliding branches: (i) mode conver sion occurs when both signs are the same, and (ii) instability occurs 
when the signs differ (e.g. ICairnslll979h . We now show that the latter situation is the relevant one in our 
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Figure 3: Dispersion relation for a cyclonic flow, with e = —0.1. The other parameters and the notation are 
the same as in Figure [2] 



problem by examining the sign of the wave momentum for KWs and IGWs in the WKB approximation. 
3.5 Wave-momentum signature 

IGWs and KWs have different leading-order approximations to their wave momentum. To see this, we 
introduce (|3.ip and (|3.3| into (|2.7|) and assume that u> is real. This gives 



M 



2N 2 m 2 
f 2 L 2 e 3 
2N 2 m 2 
f 2 L 2 e 3 



e(l - e) d\p\ 2 kco(l-s 2 + s 2 e) , 

' I .-.■}, 1 i —\ 1 1 „2.-.2\2 \P\ 



dy 



2(Co 2 - 1 + e)(l - s 2 Co 2 ) dy {Co 2 - 1 + e)(l - s 2 Co 2 ) 2 
M, (3.25) 



where the last line defines the dimensionless wave momentum M. which we will use henceforth. For IGWs, 
the first term is negligible: indeed, in the regions where p oscillates rapidly, d|p| 2 /dy = 0(1), while in the 
possible regions where p decays exponentially, d|p| 2 /dy = 0(e _1 ) only for a range of y of size e; both types 
of regions thus contribute at 0(e) to M. This leads to the leading-order approximation 

M ^ J W^W J ~^r m2 dy forIGWs - (3 ' 26) 

Given that the denominator (Co 2 — 1) cancels with the same factor in \p\ 2 (see (|3.8|l - (|3.10|l ). it is clear that 
instability involving IGWs implies that Co changes sign. It follows that there is at least one turning point in 
the channel, as announced, since the absence of turning points (£ 2 > 0) implies that |c| > 1. Assuming there 
are turning points, the sign of M. for the two types of IGWs considered in M3.3I is then 

M < for IGW+ and M > for IGW_. 
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For KWs, the two terms in (|3.25[) have a similar, 0(e), order of magnitude. Using (|3.8p . we find that 



* 2[^ 2 (±l)-l][l-s 2 w 2 (±l)] 

Using the dispersion relation for Kelvin waves, w(±l) = ^fak/r + 0(e) in our non-dimensionalisation, and 
its consequence A(±l) = r (see (|3.15[) and [|C.2|) ). this reduces to 

4 

-M±~t|^|p(±1)| 2 for KWs, 

leading to the following signs: 

M ^ for KW+ and ^ for KW_ when e ^ 0, 

differing in the anticyclonic and cyclonic cases. 

With the wave-momentum signatures just obtained, it is clear from Figures [2] and [3] that the numerous 
intersections of branches correspond to waves with oppositely signed M.. This establishes the existence 
of many modes of instability, both for anticyclonic and cyclonic shears. The main difference between the 
differently signed shears is that instabilities involving KWs only are possible only for anticyclonic shear. 

All the instabilities are associated with the interactions of modes exponentially localised on different sides 
of the channel. Therefore their interaction is exponentially weak and, as a consequence, the growth rates 
of the instabilities and range of unstable wavenumbers are exponentially small in e, as anticipated in the 
Introduction. As the asymptotic calculations of the next section show, such small growth rates are somewhat 
delicate to capture analytically. However, the interpretation in terms of interactions of waves with oppositely 
signed M. makes it possible to predict instability robustly, without detailed calculations. 



k(l - 



A(±l)[l 



l£(±l)| 5 



4 Instabilities 



4.1 KW-KW instabilities 



We start our study of the weak instabilities associated with mode interactions by deriving an estimate 
for the growth rate of the instability that aris es through the resona nce of KWs in anticyclonic shear. This 
instability has been examined in some detail by lKushner et aL ( 1998f ) and bv lYavneh et al.l ( 2001). Because it 



is the strongest instability, with physical relevance in Taylor-Couette and accretion discs (see lDubrulle et al 
120051 ). we present here a complete asymptotic derivation of the growth rate. For the KW-IGW and IGW- 
IGW instabilities considered in ^4.2lH731 we limit the derivation to the exponential behaviour of the growth 
rate as e — > 0. The method we now describe could however be applied to these instabilities as well, should a 
more accurate estimate be needed. 

To obtain the growth of the instability, we need to reconsider the dispersion relation (13.13|) in the vicinity 
of the resonance point, taking into account exponentially small terms. Let c* and r* be the values of r and 
c at resonance. By symmetry, = 0. According to (|3.16[) (with a = 1 corresponding to the anticyclonic 
shear) , 

r = r* = 1 + O(e). 

Thus, resonance occurs on an ellipse with semi-axes 1/s and 1 in the (fe,m)-plane, and the instability region 
is an exponentially small annulus around this ellipse. It is best parameterized using the polar coordinates 
(r, 8), with 

sk — r cos 8 and m — r sin 8. 

Now, take 

c = C and r = r* + R, 

where C and R are exponentially small. This can be introduced into the dispersion relation (|3.13[) ; using the 
fact that (c = 0, r = r*) satisfy (|3. 14[> . a Taylor expansion leaves only terms that are exponentially small. In 
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the coefficients of C and R in these terms, we can approximate r* by its leading-order estimate 1. Noting 
that, in this approximation, 



A(±l) w A*(±l) ■ 
we find the dispersion relation in the form 



(s 2 - cos 2 6)R ± cos 2 9(1 - s 2 )C 



s 2 sin 2 I 



cos 2 9 



s 2 sin 2 I 



where 



(R 2 -C 2 ) -4e- 2 */ £ , 



(4.1) 



My) 



'2A*(y). 

and the subscript * indicates evaluation at the resonance point. A consistent approximation of requires 
to include the 0(e) contribution to A* in the first term of the integrand. To this end, we compute the KW 
dispersion relation to O(e) in Appendix [Cl and find that = 1 + e/2 + 0(e 2 ). This leads to 



where 



with 



and 



with 



My) d y 



(4.2) 



An = 



cos 2 9(1 - y 2 ) + s 2 sin 2 



s 2 (l — cos 2 9y 2 ) 
My) cos 2 9y 2 



1/2 



My) 

1 — cos 2 9y 2 X (y)s 2 (l — cos 2 9y 2 ) An(y) 
2 cos 2 9 sin 2 9 



dy, 



(4.3) 



My) 



cos 2 % 2 



1 — cos 2 9y 2 



The second integral has to be interpreted as a Cauchy principal value at the singularities y — ±s/cos0 of 
ho(y) when these are in [—1,1]. With this result, the dispersion relation (|4.ip can be rewritten as 



C = 



where 



R 2 - a 2 e- 2 * n / e 
2s 2 sin 2 (9e-* 



1/2 



(4.4) 



s 2 — cos 2 9 

Formula (|4.4[) is the first main result of this paper. It provides the leading-order asymptotics for the 
growth rate of the KW-KW instability (after multiplication by k) as e — > and for arbitrary s < 1). It also 
makes evident the exponential smallness of the growth rate and of the instability-band width. Its validity is 
confirmed in i )4. 41 where it is compared with numerical results. 

The minimum of ^Pq, and hence the maximum growth rate, is attained for 9 = n/2, for which ^ ~ 2. 
Thus, at the crude level of exponential dependence on e, we obtain the estimate 



log Im uj ~ — , as e 
e 



0. 



(4.5) 



for the largest growth rate Imu = fclmC. Note that because 9 = tt/2 implies that k = and hence u> = 0, 
the maximum growth rate is in fact achieved for 9 slightly less than 7r/2; this does not affect the exponential 
dependence in (|4.5[) . however (see below). 
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Estimates more precise than (|4.5[) can of course be inferred from (|4.4[) . Focusing on the limit 6 — ► 7r/2, we 
note that C depends on the relationship between s and 9. A distinguished limit is found for s — 0(cos 9) < 1. 
This corresponds to the regime with s«l and 8 — k/m — 0(1), which we term the quasi-geostrophic regime, 
since it corresponds to the quasi-geostrophic scaling implying, in particular, the hydrostatic approximation 
(k/m can be recognized as the square root of the Burger number based on the wave scale). Taking the limit 
6 -> tt/2 of P~2" j) -(|3~4" )) with k = sj cos 6 fixed then yields 



■ tan 1 k and 



a 



fe-i 



\l + k\ k + 1 



The maximum of the imagi nary part of the phas e speed is then obtained for k — ► and given by Imc ~ 
2exp(— 2/e), consistent with Yavneh et alj (2001)'s equation (35). The maximum of the growth rate Imcj 
is easily seen to be attained for k — 0(e 1 ^ 2 ) and to be a factor e 1 / 2 smaller than the maximum of Imc. 
In dimensional terms, this means that the horizontal and vertical scales are both large, but have different 
orders of magnitudes, scaling like e -1 / 2 and e _1 , respectively. 



4.2 KW-IGW instabilities 

The KW-IGW instabilities occur for anticyclonic flows through the resonance of an IGW, which has one 
turning point and is localised on one side of the channel, with a KW localised on the other side. To estimate 
their growth rates, we can consider a solution consisting of a linear combination of the IGW_ given by 
(I3.19p - (|3.20[) which is oscillatory near y = — 1, and the KW+ given by (|3.12[) . (The other combination, of 
IWG+ with KW_, has the same growth rate, by symmetry.) A calculation similar to that carried out for 
KW-KW instabilities could in principle be performed to obtain the leading-order behaviour of the growth 
rate. However, this requires the derivation of the IGW dispersion relation accurate to 0(e) involving an 
inordinate amount of calculation. We shall therefore limit ourselves to the determination of the exponential 
behaviour of Imui (that is, to the determination of the constant ^>q such that loglmw ~ — VPo/e as e — » 0) 
in the instability regions, and ignore the order-one prefactor in the expression of Imw. As in the case of 
KW-KW instabilities, is determined simply from the amplitude of the colliding modes at the boundary 
where they are exponentially small, given explicitly by exp(— vfo/e). Note that controls not only the 
exponential smallness of the growth rate but also that of the width of the instability bands. 

For simplicity we restrict our analysis to the quasi-geostrophic scaling s <C 1, 8 = 0(1). For s« 1 and 
(7 = 1, the phase speeds of colliding KW+ and IGW_ branches given in (|3 . 16[) and (|3 . 23[) reduce at leading 
order to 



1 



1 



and 



-1 



1 

fc 2 



1 



1/2 



respectively. The corresponding resonance condition 



1 



1 



-+ To+-7 =2, 



that is, 



If TO 



2 V 771 — 1 



TO- 



1/2 



1/2 



with to > 1 



(4.6) 



defines a curve in the (k, to) plane in the vicinity of which instabilities are concentrated. For KW-IGW 
instabilities, since there is a single turning point y_ in the channel, is given as 



(4.7) 



*o= / Hy)dy. 
Jy- 

The integrand X(y), given in (|3.7p . can be approximated by 

A(y) = km [(y + - y)(y - y~)] 1/2 , 



(4.8) 
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1 



1 



1/2 



with y± reducing to 

y± = c± 

Introducing (|4.6|) and (|4.8|) - (|4.9|) into (14. 7[) gives the expression 

1 



3 - 2/m 
-1 



(4.9) 



*n = 



(2m- 1 



f ( 7T + 2 sin" 1 — - — ^ + 4(m(m - 1)) 1/2 
\ 2m — 1 / 



8[m(m-l)]Va 

The maximum growth rate of the KW-IGW instability is given by the minimum value of \1/oj found to be 



*o = 2.80 ■ 
Thus we obtain the asymptotics 



log Im lu 



for fc = 1.04--- and m=1.30- 
2.80 



as e — > 0, 



(4.10) 



(4.11) 



for the growth rate of KW-IGW instabilities. Comparison with (|4. 5[) then indicates that these are consider- 
ably weaker than the KW-KW instabilities. 



4.3 IGW-IGW instabilities 

We now consider the instabilities that result from the resonance between IGWs. These are particularly 
important for cyclonic flows since they provide the only mode of instability in this case. In fact, as can be 
expected from the leading-order dispersion relation (|3.23p . the dominant behaviour of these instabilities is 
unaffected by rotation, so that the exponential dependence on 1/e is identical for anticyclonic and cyclonic 
shears. What differs between the two cases, however, is the order-one prefactor which we do not estimate 
analytically. 

IGW-IGW instabilities occur when a solution p_ of the form (|3.19j) - (|3.20p is resonant with its counterpart 
p + . The modes have then two turning points y± in the channel, leading to the necessary condition r > 
(1 + 5 2 ) 1 / 2 for the instability. We now estimate the factor vf^ controlling the exponential smallness of the 
instability growth rates. As in the previous section, we restrict our attention to the quasi-gesotrophic scaling 
s < 1 and S = 0(1). We furthermore consider only the strongest IGW, associated with the (symmetric) 
resonance of the gravest (n = 1) IGW modes, and for which c = to all orders in e. The resonance condition 
is therefore 

1 1 

Since for n = 0(1), the two turning points are y± = ±1 at leading order in e, 4"o is computed as 

q ' o = km J y i(y+ - y)(y - y-)] 1/2 d v = -p- 

The minimum value is therefore 

* = for k — m — %/2, (4.12) 
and the exponential scaling of the growth rate given by 

7T 

loglmw ~ , as e — > 0, (4-13) 

for both anticyclonic and cyclonic flows. This is exponentially smaller than the growth rate for either the 
KW-KW or the KW-IGW instabilities g3|) or (jiTTjl . 
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Figure 4: Imaginary part of the phase speed Im c as a function of r (in linear-logarithmic scale) for the 
KW-KW instability in anticyclonic flows with s = 0.1. Numerical (dots) and asymptotic (solid curves) 
results are compared for e = 0.3, 0.4, 0.5 (with increasing Imc), and 9 = n/2 (i.e. k = and m = r). 



4.4 Numerical computation of growth rates 

We now present comparisons of the growth rate, or rather Imc, computed numerically with the asymptotic 
results of § §4.1H4~3l The numerical method employed is that described in §3.41 where Re c was considered. 
For the small values of e examined here, Im c is very small and the bands of unstable wavenumbers are very 
narrow, so that very fine resolution in y is needed to capture Im c accurately. In order to ensure high accuracy, 
we successively double the grid resolution until results are unchanged to at least four significant digits. This 
required grids of sizes ranging from about 250 mesh points for strong or moderate instabilities, to as many 
as 16 000 mesh points for very weak instabilities. This may be improved upon by using nonuniform grids 
with high resolution only in regions where the solution changes fast. The search for the bands of instabilities 
in (k, to) is quite delicate, but made possible by the excelllent approximations afforded by the asymptotic 
results. 

We start by considering the KW-KW instability of anticyclonic flows. Figure [4] shows Imc as a function 
of r for 9 = n/2 and e = 0.3,0.4,0.5 in the instability bands. The dots represent numerically computed 
values; the solid line are computed analytically using (|4.4|) . Note that we only know r* to algebraic accuracy, 
while the bands are exponentially narrow. Hence, we use the numerical results for determining r* — the value 
of r for which Imc is maximized. The narrowing of the instability band is clearly exhibited in the figure, 
and the small-e analytical approximation quickly converges to the numerical results as e becomes small. The 
dependence of Im c on 9 is illustrated by Figure [5] which compares numerical and asymptotic estimates for 
the maximum value of Imc as a function of 9 for s = 0.1 and e = 0.3, 0.4 and 0.5. The value of Imc in the 
quasi-geostrophic scaling s -C 1, S = O(l), that is, the limit 9 — ► n/2, is also indicated. The Figure confirms 
the accuracy of the asymptotic estimate and shows the rapid decrease of Imc as 9 decrases from n/2. 

Our results for all the types of instabilities are summarized by Figure [6l This compares asymptotic 
and numerically computed values of Imc as a function of l/\e\ for KW-KW, KW-IGW and IGW-IGW in 
anticyclonic flows, and IGW-IGW instabilities in cyclonic flows. The values of Imc displayed correspond 
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Figure 5: Maximum of Imc as a function of 9 (in linear-logarithmic scale) for the KW-KW instability of 
anticyclonic flows with s = 0.1. Numerical (dash-dotted curves) and asymptotic (solid curves) results are 
compared for e = 0.3,0.4,0.5. The limits of Imc as 9 — > 0, corresponding to the quasi-geostrophic scaling 
s <C 1 and S = 0(1) are also indicated. 
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Figure 6: Maximum of Imc as a function of l/|e| (in linear- logarithmic coordinates) for all the instability 
mechanisms examined in this paper, both for anticyclonic (e > 0) and cyclonic (e < 0) flows. The asymptotic 
estimates (solid lines) are compared with numerical results (symbols) for s = 0.1. 



to the maximum over m and k for fixed s = 0.1. For KW-KW instabilities, the asymptotic estimates are 
obtained from (|4~4|) . For KW-GW and IGW-IGW instabilities, we use (|4~TT1) and (|4TT3]) . respectively. These 
give Im c only up to a multiplicative constant which we fix by matching the asymptotic and numerical results 
for the smallest values of |e| shown in the figure. In the linear-logarithmic coordinates used, the numerical 
points line up with the predicted straight lines for larger |e|, thus confirming the validity of the asymptotic 
analysis. Further support is provided by the fact that the values of k and m for which Imc is maximised 
are close to the estimates (14. 10|) and (|4.12[) . Evidently, the match between the numerical and analytical 
results is quite good even for e moderately small. We see that the instabilities become substantial for e w 1, 
especially KW-KW instabilities. Observe that, as predicted by the analysis, the decay of the growth rate in 
IGW-IGW instability as e becomes small is the same for cyclonic and anticyclonic flows, and yet the growth 
rates of cyclonic flow are smaller by a factor of about 20. Thus the O(l) prefactor in the asymptotics of 
Imc for IGW-IGW instabilities, ignored in (|4. 1 3[> . turns out to be numerically very different for anticyclonic 
and cyclonic flows. The smallness of this prefactor in the cyclonic case means that the instability remains 
exceedingly weak even for ew 1, and likely irrelevant in many physical situations. 

5 Discussion 

This paper examines the linear stability of a horizontal Couette flow of a rapidly rotating, strongly strati- 
fied, inviscid fluid. The main conclusion is that the flow is unconditionally unstable: unbalanced instabili- 
ties, associated with linear resonances between Kelvin and inertia-gravity waves, occur for arbitrarily small 
Rossby numbers e = A//. The growing perturbations have small horizontal and vertical scales, with typical 
wavenumbers or spatial-decay rates of the order of e _1 . Physically, it is easy to understand why asymp- 
totically small scales are a key ingredient of the instabilities. The phase locking between different waves 
which underlies the instability mechanisms requires the wave phase speed to be comparable to the basic flow 
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velocity, and this only occurs for small-scale waves. The need for small vertical scales also explains why 
the instabilities examined in this paper have no direct counterparts in shallow-water flows; these are stable 
for small enough \e\ because of the inherent limitation in vertical structure imposed by the shallow- water 
approximation. 

Our conclusion that the rotating stratified Couette flow is always unstable is of course in sharp contrast 
with the one that may be drawn from balanced models. Regardless of their accuracy, which can be any 
power e™, they predict the stability of flows without inflection points such as the Couette flow. There is no 
contradiction, however, since the growth rates found for the unbalanced instabilities are exponentially small 
in |e|. In practice, this exponential dependence means that the instabilities are exceedingly weak when |e| is 
small, but can become important rather suddenly as |e| increases towards 1 and beyond. If the instabilities 
are to play a significant role in the breakdown of balance in geophysical flows, this will therefore be in a 
manner that is extremely sensitive to the Rossby number. 

In the literature, most attention has been paid to anticyclonic flows, and in particular to the coupled 
Kelvin-wave instability occuring in these flows. Our results clarify that cyclonic flows are also unstable, 
through an instability mechanism involving coupled inertia-gravity waves. This mechanism is also active 
in anticyclonic flows where, along with the instability mode mixing Kelvin and inert ia-gravity waves, it 
provides an alternative to the well studied instability due to Kelvin- wave resonance (see lYavneh et al. 2001 



iMolemaker et all 2001 ). The focus on anticyclonic flows and Kelvin- wave instabilities is justified in practice 
by the fact that the associated growth rate is much larger than those of the other instability mechanisms, 
exponentially larger in fact in the limit e — > 0. The instability of the cyclonic flows is especially weak. This 
weakness is not completely accounted for by the exponential dependence on 1/e, since this is the same for 
both anticyclonic and cyclonic flows whilst the growth rates obtained numerically are very different. We 
conclude, then, that the exponential dependence and the O(l) prefactor conspire to make the instability of 
cyclonic flows extremely weak, even for moderate |e|. 

The WKB approach used in this paper could be extended to examine the instability in more general 



rotating stratified shear flow s. Obvious applications are the stratified Taylor-Couette flow (jYavneh et al 



2001 . Molemaker et al.ll2001 ). which differs from the prob l em studied here by th e presence of curvature terms, 



and the stability of accretion discs ( Riidiger et al. 120021 iDubrulle et al. 1 120051 ). Additional physical effects 
that it would be of interest to study include different boundary conditions (in particular the case of infinite 
domains for which no Kelvin waves exist), viscous and thermal dampin g, and non-zero potential- vorticity 
gradients, leading to the existence of critical levels for neutral modes (cf iBalmforth 19961 ). 
JV was funded by a NERC Advanced Research Fellowship. 



A Conservation laws 

Let 

M = {ud z p-wd x p)/N 2 . 
Denoting integration over the periodic domain in x and z by 



(•)=//■ dxdz > 

we compute 

N 2 d t (M) = (d t ud z p~d z udtp-d t wd x p + d x wdtp) 

= ((/ — A)vd z p — d x pd z p — N 2 d z uw + d z pd x p) (A.l) 
= -N 2 d y {uv), 

where we have used integration by parts and periodicity extensively, and, for the last line, q — and the 
incompressibility equation. The conservation for the quadratic wave momentum (or pseudomomentum) 

M = [ [ I (ud z p - wd x p) Axdydz/N 2 
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follows by integration in y, using the boundary conditions v = 0. 

The perturbation energy £', with density |u| 2 /2 + p 2 /(2N 2 ), is not conserved but satisfies 




Airy dxdydz. 



Integrating by parts the right-hand side and using (|A.1[) gives a conservation law for the wave energy (or 
pseudoenergy) 

Note that the conservation of both Ai and £ can also be derived from the exact conservation laws for 
momentum, energy and potential vorticity for the full system, that is, basic flow plus perturbation. 



B Equation for u 

In §2, the eigenvalue problem satisfied by normal-mode solutions is for mulated as the second -order differential 
equation (|3.4|1 for p and i ts associated bounda ry condition ()3.5|1 (cf. iKushner et al. 1998f ). An alternative 
formulation, employed by lYavneh et al.l (|200lh . uses u instead of p as the dependent variables. It has the 
advantage that the removable singularities that appear in (|3.4p are absent. For completeness, we record this 
alternative formulation as 

2 / l-.s 2 cI> 2 V / k 2 (l-s 2 Cu 2 )+m 2 (l-e-LU 2 ) 2e(l - e)s 2 k 2 m 2 u, 2 \ . _ 
6 V K U ) \ K + K 2 ) U ~ ' ^ A ' 

where 

K = (1 - s 2 uj 2 )k 2 + (1 - e) 2 m 2 . 
The associated boundary conditions are 

fc(l — s z u) z ) 

This is the formulation used for the numerical computation of the normal modes. 



C Kelvin- wave dispersion relation 

In this Appendix, we derive the dispersion relation for KWs accurate to O(e), as is necessary to obtain the 
leading-order asymptotics of the KW-instability growth rate. 

The dispersion relation for KW± valid to all orders in e are given in (|3.14|) . It is solved at leading order 
in i j3.2l to give (|3.16|) . At the next order, we find the two equations 

± a Cl A(±l). 9± (±l) + c (±l)<4(±l) = 0, (C.l) 

which allow the determination of the 0(e) contribtion to the frequency lo\. Note that the contributions of 
the O(e) terms neglected in (|3.11|l - (|3.12|) cancel in these two equations when (|3.15[) is taken into account. 
Equation (|3.9[) can be used to express the derivatives of g± ; the following results are therefore useful: 

(C.2) 



A(±l) 


= r, 


A'(±l) 


±a(l - s 2 


2A 


m? 


*wo(±l) 




1-^(±1) 


r 2 — k 2 ' 


Tcrh{±l) 
2A(±1) 


= qForQ- 



" 2 — k 2 



{l-s 2 )k 2 r 2 

Ci = . 
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Using these and (|3.9[) , (jC.ip gives the first-order correction to the frequencies (|3.15[) , 
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